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ABSTRACT 

The behavior of weak electrostatic waves in a collisionless mag- 
netoplasma supporting a steady large amplitude whistler wave has 
been studied. All waves are assumed to propagate parallel to a uni- 
form background magnetic field, B 0 . In the presence of the whistler 
wave fields each particle executes an oscillatory motion parallel to B 0 , 
in addition to a translation along B 0 and transverse motions. This os- 
cillation causes the Landau resonance to be replaced by a series of 
new resonances between particles and the electrostatic modes. A dis- 
tribution function for the perturbed plasma is constructed by solving 
the Vlasov equation, linearized in the electrostatic wave amplitudes. 

A dispersion relation is obtained and solved approximately for the 
growth/damping rate of the perturbations. Growing electrostatic 
modes are found to be approximately uncoupled. Trapped particles 
have a strong influence on the stability of the system. 
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STABILITY OF A STEADY, LARGE 
AMPLITUDE WHISTLER WAVE 

I INTRODUCTION 

The existence of steady, arbitrary amplitude whistler mode solutions to 
the Maxwell-Vlasov equations was demonstrated in 1966 by Lutomirsld and 
Sudan 1 . These authors considered right circularly polarized waves propagating 
along a uniform background magnetic field, B 0 , in an infinite collisionless 
plasma. It was shown that no dispersion relation exists for such waves. Un- 
damped, finite amplitude whistlers with arbitrarily chosen frequency co 0 and 
wave number k 0 are therefore possible, in principle. 

In order to decide which of these modes are likely to be observed in real 
plasmas, one must be able to answer two questions. First, does a physical 
mechanism exist which "constructs" the mode in a plasma? Second, is the re- 
sulting equilibrium stable enough to persist for a reasonable time? In this 
paper we consider the question of stability. Specifically, we investigate the 
stability against the growth of electrostatic perturbing waves of a plasma sup- 
porting a large amplitude steady whistler. Of particular interest is an electro- 
static instability associated with the distribution of trapped particles. 

There is a second motivation for studying the stability of undamped whistlers. 
A recent theoretical paper 2 describes a process which produces, in the absence 
of instabilities, an approximate whistler mode equilibrium. In the case consid- 
ered in Reference 2, a right circularly polarized transverse wave propagates 
along B 0 in an infinite collisionless plasma initially in thermal equilibrium. 

The wave damps exponentially for a short time after it is launched at t = 0. 
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The resonance damping mechanism acting on the wave rapidly saturates, due 
to the phase mixing of resonant particles. The wave amplitude then executes 
damped oscillations about a finite value, and the system approaches a steady 
state. The presence of instabilities may modify this process. One can avoid 
the difficulties inherent in a stability analysis of the complete, time dependent 
system by considering just the asymptotic steady state. This can be done in a 
straightforward manner, and yields some insight into the nature of the instabil- 
ity mechanisms involved in the dynamic process. 

We begin in Sec. II by linearizing the Vlasov equation for the electrons in 
the perturbation amplitude. Ion motions are neglected and all waves are as- 
sumed to propagate along B 0 . The solution to the linearized Vlasov equation is 
constructed in terms of unperturbed single particle trajectories. (The equations 
of motion for a single particle in B 0 and the whistler fields have been solved by 
Roberts and Buchsbaum , 3 and others. 1 ’ 2,4 ) The perturbed distribution function 
is used to obtain a dispersion relation for the Fourier components of the per- 
turbation. 

It is known 1 ' 4 that the fields of a finite amplitude whistler wave trap par- 
ticles, and that the motion of both trapped and untrapped charged particles in 
these fields consists of an oscillation parallel to B 0 , a translation along B 0 , and 
transverse motions. One finds that the parallel oscillation part of the unper- 
turbed particle motion, which is greatly enhanced for trapped particles, leads 
to new resonances with the electrostatic perturbing waves. 

One also finds that the perturbing waves are approximately uncoupled, due 
to the symmetry of the unperturbed wave-plasma system. If we assume that 
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the whistler phase velocity a> 0 / k 0 is less than c and work in the wave frame, 
there is no electric field and the magnetic field lines are helices. 

B(z) = B t (e x sin k Q z + e y cosk 0 z) + B fl e z 

The unperturbed electron distribution function also has helical symmetry, i.e., 
the electron distribution at any point z t differs from that at another point z 2 only 
by a rotation in velocity space through an angle i p - k Q (z 2 - Zj). The electron 
density is therefore uniform throughout space. This is in contrast to the situ- 
ation encountered by Kruer, Dawson, and Sudan 5 , and Goldman 6 , in their study 
of the stability of a large amplitude electrostatic wave. In that case there is a 
strong periodic spatial variation in the unperturbed electron density, and this 
leads to the coupling of an infinite number of discrete perturbing waves. 

Finally, in Sec. Ill, an approximate expression for the growth/damping rate 
is derived and studied, with particular emphasis on the effects of trapped 
particles. 


II THE DISPERSION RELATION 

The wave frame electron distribution function for a magnetoplasma sup- 
porting a large amplitude, parallel propagating steady whistler is a function of 
the single particle constants of motion 1 ' 4 v and W, where 

v = ( v i| 2 + Vj 2 ) 54 (la) 

1 / fi oV ^1 

W " 2rl + kJ + V 1 sin ( k o z + s^) ( lb ) 


The quantities v |( , v L , and yp describe the wave frame velocity 
V = V L cos I pe x + V x sin I pe y + v |( e z 
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and 


0. 0 = | e | B 0 /m, fi, = |e|B,/m. 

A small electrostatic perturbation of the form E|| (z, t)e 2 produces a cor- 
rection to the distribution function (f - f + f t ) which may be computed by solving 
the linearized Vlasov equation. 

3 3 e e 3 f (v, W) 

— * V || 5 — v x B(z) • V f. - -- Eii(z, t) — k" — 

' t 11 3 z m v ' v J l m 11 v ’ ' 3 v || 

( 2 ) 


e _ 

v ll 3 f / 


*\j 

\ 3 f 

E|i 
m 11 

<1 

<1 

■ 4 - 


1 3W 


This equation may be solved formally by introducing the function z(z 0 ,v 0 , t), 
which describes the zero order motion along the z axis of a particle initially 
at z 0 , v 0 . The functions z 0 (z, v , t) and v 0 (z, v , t) are, like v and W, constants 
of the unperturbed motion and solutions of the homogeneous equation associated 
with (2). Thus 


f, 


^lo( Z 0’ V o) m || j' E||( z (t'). t') 


( 2 a) 


I d »(t') 3 f.(v, W) [ dz ( t ') % \ 3 f ( v, \V) 

vdt 1 3v ydt' k Q J 3 W 


dt' 


If one considers only growing disturbances, 

f E,|(z(t'), t')[ " ]df 

J -CD 


e f 0 
m 


( 2 b) 
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First order corrections to the perpendicular current field and the charge 


density can be obtained directly from f 

J i, = " n e e///v 1 fjd 3 v p, -• - n e ej'/jf t d 3 v 


The current density J tJ , which is in general nonzero, leads to time varying 
magnetic and induced transverse electric fields. The omission of these fields 
from (2) is justified when the perpendicular thermal velocity V T is much less 
than c. To see this, we estimate that J, l ~ V T pj ~ e 0 V T k E M , where k is a 
typical wavenumber in the spectrum of Eh (z, t). After using Maxwell's equa- 
tions to relate the transverse field corrections E lL and Bj^to J^, and introducing 
the characteristic frequency co i and wavenumber k t for E^, B^, we find that 


* \ k V T V 

— L ^ j— < < i 

k l c 2 + aif/k* 


Z±L 
E .I 


k_ V T («iA|) 

2 + a)j 2 /k 2 


<< 1 


l c 


The charge density p y produces an electric field which must be identical 
with E H (z, t), for self consistency. If we express E| ( (z, t) in the form 

E,| (z, t) = f E ll( k) e‘<— t) dk, 

J -CO 

construct the charge density by integrating (2) over velocity space, and apply 
Poisson's equation, we obtain 
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where on(k) = &»(k) - k w Q /k 0 and o> p 2 = n e e 2 /<- 0 m. a>(k) is the complex frequency 
of E llk (t) in the lab frame. 


As noted previously, the unpei'turbed motion of a particle along the z axis 
consists of a pure translation plus an oscillatory motion. If t + (z 0 , v Q ) is any 
time at which the particle reaches the upper turning point of its oscillatory 
motion, the function z (z 0 , v 0 , t) may be written in the form 

z = z 0 t v,(v,W)t + [S(v, W, t - t + ) - S 0 ] (5) 


where S(v, W, t - t + ) is periodic in time with period T (v, W), S 0 is S (v, W, - t + ), 
and Vn (v, W) is dz/dt averaged over a period. The functions Vu , S, T, t + are 
known, in the sense that they can be constructed in a straightforward way from 
the exact, closed form expressions for dz/dt obtained by Roberts and 
Buchsbaum (Ref. 3), or from the simpler approximate expressions given in 
Ref. 2. The time dependence of exp (i k z (t)] , which appears in (4), can be 
conveniently exhibited by expanding exp [ i k S (t) ] in a Fourier series, so that 


exp 


kz(t)|= exp ik(z Q + V || t - S Q ) j ^ A n exp J^27rin(t - t + )/T 


(6a) 


where 


A n (V,W,k > = f(v!\V) GXP 

•'n 


T r (t - t + ) 

ikS(t - t + ) - 2win ^ 


d(t - t + ) (6b) 
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After noting that [exp(i k z)] dz/dt is (i k)" 1 d[exp (i k z)]/dt and ap- 


plying (6), one is able to write (4) in the form 


o = r E,<k)< 


e i (ki-w't) |< 


-1 (j', t 


' JOT d 3 v 


5 ik( *o“ s O 


J -oo n 


t “ 


exp 


i k v,| t + 7m\ 


d f > dk 


where 


I„(v. W,k) 


°0 3? / 2rm\ / 13? 3 f\ 
k 0 3W + \ V|1 + kT / 3v + 3 \vj 


( 7 ) 


After integrating over t\ using (5) to express (z 0 - S 0 ) in terms of current 
dynamical variables, and expanding exp (- i k S), we obtain 

0=1” E„(k) o 1 *> | k 


(A v 

k 




exp [ 2v i (n - m) (t - t + ) / T] 1 


(v., + 27rn/kT - oj< /k 


- f dk 


I 


The terms in the double summation above for which (n - m) is not zero phase 
mix away for large t, when integrated over velocities. Note also that the quantity 
in curly brackets in (8) is independent of z. The implicit z dependence in 
W (v,| , v.l, k 0 z+ 4>) is removed by the integration over (d 3 v - v A dv x dv h d> p). 

It follows that the Fourier components Eh (k) are uncoupled, and the disper- 
sion relation is 
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0 


( 9 ) 


a;,: 


k2 


mZ k< v - w >i 2 


I„(v,W,k) cl 3 v 

(V|| + 2?/n/k T - w/k) 


This expression was derived for growing modes, hence it is understood that the 
integrals are to be taken along the real velocity axes when the frequency has 
a positive imaginary part. The integration contours are deformed in the usual 
maimer to construct the analytic continuation of (9) in the lower half w-plane. 

It is easily verified that (9) reduces to the dispersion relation obtained by 
Landau in the limit 


Q. - 0 (A n - S Mi0 ; v n - V|| ; I Q - 3f/3v ). 

Ill ESTIMATION OF Im M 

There is an alternative form of (9) which provides a convenient starting 
point for a study of the resonant amplification or damping of perturbing waves . 
After a little algebraic manipulation, and after using the fact that 2 1 A n | 2 = 1 
[ from (6) ) , (9', may be written in the form 


e (k,w) 



i a f ii 
v 3 v *" 3 W + 



|A n | 2 I(v,W, w/k) 

( V || + 27rn/kT- w/k) 


d 3 v = 0 (10a) 


where 


I (v. W, «•/ k) 


o' 3 f oji 
k v 3 v ] k H k Q 


3_f 
3 W 


(10b) 


When the growth/damping rate is small, the imaginary part of (10a) may be 
expanded in powers of / co p , where co i = Im ai>. Alter neglecting terms 

higher than first order and solving for w l , one obtains 
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) 


0 t <1 (V, W, «./ k)> k /0(k,6J r ) 


where w r is the real part of the frequency of the perturbation in the lab frame, 


K \ s n 


ST' 

— J/J( )^|A n (v,W)| 2 5(Un)d 3 v 


Q(k, co r ) = o; r k — 6 r (k. <a r ) 


^ 2 Drrr Jl3f 9? yl A n| 2l l,3 /11H V 

ir( k,» r) = z,— |d 3 v (lid) 


u n = v H ( v, W) + 2wn/kT - <u'/k 


Similarly, the real part of (10a) yields e r (k, w r ) % 0 when u i is small, and 
this determines the function « r (k). The quantity Q (k, co r ) determines the sign 
of the energy associated with the wave E k in the lab frame . In a Maxwellian 
plasma Q is positive for all values of k, but negative energy waves are possible 
in media sufficiently far from thermal equilibrium ("anomalous dispersion") 7 . 

The growth or damping rate of perturbations on a given whistler equilibrium 
can, in principle, be determined by direct application of (11), but the computa- 
tions involved are clearly formidable. It is therefore helpful to obtain some 
qualitative or semiquantitative information about factors influencing the stability 
of the system. To simplify the problem, we assume that B t is small compared 
to B 0 and study effects of lowest order in e, where 

e = b/Bq = n,/n 0 . 




(’"WSi, 








The Cyclotron Resonance 


We begin by considering the unperturbed particle orbits. The unperturbed 
path of a particle in the two dimensional phase space with coordinates k 0 z + , 

v || can be determined by solving (lb) for k Q z + 0 and plotting k Q z + i p versus v,, . 
A typical set of such paths, for particles having the same v and various values 
of W, is illustrated in Fig. 1. Note that the angle 90° + k Q z + < Jj, which is the 
angle between and v L , oscillates about zero for particles following the closed 
orbits. These particles are said to be "trapped" by the whistler. Untrapped 
particles close to the separatrix and trapped particles interact strongly with the 
wave and are termed "cyclotron resonant." (w 0 " k 0 V, LAB FRAME ~ ~ ^ 0 V l> ~ 
n 0 for cyclotron resonance.) 

The following results are based on the approximate unperturbed particle 
orbits developed in |lef. 2. For trapped particles, 

l/n 0 T(v, W) = e' A [(k 0 v/n 0 ) 2 - l] ' V 4 K ( 1/r ) + 0(e) 


v„(v, W) = -Q 0 /k 0 +O(e) 


(12a,b,c) 


S(v, W, t + ) t) - p 


8 Vi q"i7r) sin[27r(2m + 1) (t - t + )/T] 


rri= 0 


(2m f 1) (1 + 


+ 0(e*) 


Here K (1/r) is the complete elliptic integral of the first kind, and 


q ( 1/r ) - exp { - -v K [ ( 1 - 1/r 2 /'*]/ K ( 1/r)} 


r (v, W) = ( 2e) 1/ 


VA 2 


- i 


a. 


W + e 


'k 0 v \ 2 


a 


- i 


L 


K I 


(13) 
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For untrapped cyclotron resonant particles, 

l/0. 0 T(v, W) = e* [(k 0 v/Q 0 ) 2 - 1 ] 14 / 2r K (r) + 0(e) 

V|| (v, W) = - O 0 /k 0 + 2wo-/k 0 T(v, W) + 0(e) (14a, b,c) 

<° qm 

S(v, W, t t , t) = -jp / — <0 — sin [27rm(t - t t )/T] * 0(e“) 

•fei "I 11 Vo) 

where a = sgn(V|| + n 0 /k 0 ). Finally, for non- cyclotron-resonant particles, 

l/n 0 T(v, W) = e* [(k 0 v/n o ) 2 - l] ,/4 Ar + 0(e) 

v,| (v, W) = - ft 0 /k 0 + 27j-cr/k 0 T(v, W) + 0(e) (15a,b,c,d) 

S(v, W, t + , t) = 0(e) I A n (v , W) | 2 = S ni0 + 0(e) 

The parameter r (v, W) "classifies" the unperturbed orbits according to type, 
r (v, W) is greater than one for trapped particles, with r - 00 for particles in 
the center of the nested set of closed orbits in Fig. 1. r (v, W) has the value 
unity for particles on the separatrix and is less than one for untrapped particles, 
with r tz O (e' A ) for non-cyclotron resonant particles. 

Locus of Electrostatic Resonances 

Resonance between particles and a perturbing electrostatic wave occurs 
when the Doppler shifted frequency of the perturbing wave, as measured in a 
frame moving with the average velocity v„ of a particle, matches a harmonic 
of the particle's oscillatory motion. It follows that 
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Y(v,r,k) = [k V|| (v, r) - coj (k)] T(v, r)/27f 


n 


( 16 ) 


at electrostatic resonance, where n = 0, ±1, ±2, . . . . With the aid of (12) — (15), 
we depict in Fig. 2 the variation of [kvn - toj ] T/2 v with 1/r (v, W), for a fixed 
value of v and for particles above and below cyclotron resonance velocity. The 
intersections of the curve Y (v, r, k) with the horizontal lines Y = 0, ±1, ±2, . . . 
mark the locations in phase space of electrostatic resonances. Note that the 
higher harmonic resonances tend to cluster in the vicinity of the cyclotron 
resonance. This is due to the singularity in T (v, r) at r = 1. (T (v, r) ~ 
r K (r) rln [ (1 - r 2 ) ,/S ] for untrapped particles, etc.) 

It follows from (12a), (12b), and (16) that particles on the flat portion of the 
curve near <r/r ~ 0 in Fig. 2 resonate with the perturbing wave when 

(v 2 - n 0 2 /k 0 2 )/v T 2 = |( W . +ka 0 /k 0 )/n(k 0 n 1 v T )«| 4 . 

These particles oscillate coherently, which greatly enhances the effect of the 
resonance. The frequency (kgO.j V T ) ' A is the "bounce frequency" co B for typical 
trapped particles at the bottom of the magnetic wells, i.e., 2rr/T (v, r) -* « B 
when r - v 2 - Q 0 2 /k 2 % v x 2 - V T 2 . 

Trapped Particle E ffects 

In view of the preceding, one expects cyclotron resonant particles, es- 
pecially trapped particles, to play an important role in the interaction between 
the whistler and the perturbing waves. The effects of trapped particles are 
most prominent when the Doppler shifted frequency of the perturbing wave, as 
measured in the frame moving with the average velocity v u = - fi Q /k 0 of the 
trapped particles, is of the order of the bounce frequency co B = (k 0 Q j V T )' A . 
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It is therefore convenient to introduce the frequency w", where &j" (k) = 

(k) + kfi 0 /k 0 , and to consider k chosen so that \a>" (k)l 'V. O (w 0 ), 

It is instructive to solve (10b) and (11) for < 3 f/3W > k . We assume that 

♦ 

the perturbing waves have positive energy in the lab frame (Q > 0) , and that 
co 0 < fi 0 . The latter assumption is reasonable in the small e limit, where one 
expects the whistler dispersion relation obtained from linear theory to hold 
approximately, but is not necessary in general. One finds that the system is 
unstable (w. >0) if 



or 




(17a, b) 



< (Q„ ~ "o)/ k o 


The quantity in square brackets in each of these inequalities is positive. 

We first consider the case in which < 3 f/v3v > k is positive for some 
choice of k. If one considers trapped particles only, it is possible to show (Ap- 

r\j 

pendix A) that < 3 f/3 W > k and < 3 f/v3v > k are even functions of w", when 
I co" | ~ O ( w B ) and to lowest order in 6 . Given this symmetry it is easy to con- 

'V/ 

vince oneself that either (17a) or (17b) is satisfied if < 3 f/v3v > k is positive, 
indicating that trapped particles are destabilizing. 

The case in which < 3 f/v3v > k is negative may be more important physi- 
cally, since (3f/v3v) w ~ (3f/v x 3 v x ) v |( in the small e limit. In this case (17) 
implies that the plasma is unstable if both of the following conditions are 
satisfied: 
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and 



(18b) 


The fact that < d[/dVI > k is even in guarantees that (18b) will be satisfied 
for some choice of k. A semi -quantitative trapped particle instability condition 
may be obtained from (18a). We estimate (3 f/v3v) w ~ - f/V T 2 and Bf/3W ~ 
/Sf/AW. AW is the change in W in going from the separatrix to the bottom of a 
well. For typical values of v this is Qj, V T /k 0 . /3 scales the fractional varia- 
tion of f over the trapping region, in the direction normal to the surfaces of 
constant W in phase space (Fig. 1). After setting I | ~ (k 0 fij V T ) 4 , one finds 
from (18a) that trapped particles are destabilizing if 


l/3| 


> 

'V, 


k (%_ 

k 0 \ k 0^T 



Special Cases 


For example, consider a distribution in which the trapped particles are 
strongly concentrated near the centers of the nested sets of closed orbits (one 
set for each value of v), and approximate the distribution function in this region 
of phase space by 


(2n 


eT' 


/n„V T 2 ) 


-v_2/V.2 


S ( V II f V k o) S(i// + k 0 z + 77/2) (19a) 
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where n cT is the density of trapped particles. It is shown in Appendix B that 
this is equivalent to 





k o‘Vr H ( v2 ~ Q o /k o) 

2 ' /TQ 0 n e 7 ?( V ) V T 


exp {- 


(v 2 - »0 2 /t<0 2 ) 

v T 2 



(19b) 


where H (x) is the unit step function, and 


V (v) 11 e 1 ' 4 


k 0 vN2 




e- 


lim sf— — — 

r 


R^ 00 \r2 p2 


Substitution of (19b) into (10b) yields, to lowest order in e, 


"" n eT k 0 3H ( v2 “ W) 

I £ exp 

477 k n e V T 2 r) 3 Q 0 3 


(v 2 - % 2 /K 2 ) 


V2 

V T 




o (20) 


where 


S' ( 1/r 2 - 0) a B 8/3 r’2, 


It follows from (16) and (12b) that the n = 0 resonance condition for trapped 
particles is -<tj"/1c = 0. Since I (v, W, w') vanishes when w" = 0, the n = 0 
resonance makes no contribution in this case. For very large r q (1/r) ~ (4r)~ 2 . 
After using (12c) in (6b) and retaining only terms of lowest order in 1/r, one has 

l A ±i ( v > w )| 2 ~ C 1 / 1- ) 2 ( k / k o ) 2 r ~ oo ( 21 ) 

The higher harmonic oscillator strengths go to zero faster than (1/r) 2 , and 
these terms do not contribute to <I> in this case. 

Equations (20) and (21) may be used in (11) to obtain a growth rate for the 
perturbations. The integrations are most conveniently performed in terms of 


15 


/ 


the variables 8 r, v, and £ = ^ 7T+ ^o z + '/'j • ° ne obtains from (la,b), (lie), 

(12a, b) and (13) the results 


where 


8k K ( 1/r ) (v 2 - n 2 /k 2 )* 

8( "" ) = .mKk.n,). 8<v ’ ‘ v ' ?) 8 "">-<“"> (22a) 


Q 2 

0 u o 
■ 2 = + 

n k 2 
K 0 


4 K (1/r) a” 


2ttx\ 


Ct B 


V 2 

T 


and 


d 3 < 


cr(x) = sgn x 

4O 0 i7(v) (1/r 2 - sin 2 ^)"Vk 0 r 3 


vdv dr df 


(22b) 


The variable £ appears only in the expression for the volume element d 3 v, and 
leads to an integral of the form 

/*+ s in - 1 ( 1/r ) 

| (1/r 2 - sin 2 £)-K d£ = 2 K ( 1/r ) 

s i n ~ 1 ( 1/r) 

in (lib). The remaining double integral contains two S functions and is easily 
evaluated. If the approximation 


v ,r 


“e 


3 k ’Vf 

e ,«vk) * i-- 2 1* 


is used to evaluate Q ( w r , k) , one obtains for co i 


co. 

J ro 


7TO) r ^o! T w"/a)gj 2 e /Wb) sgn(w") 


(23) 




where ^ T 2 = n eT e 2 /s m ' The factor exp {- (w"/c<j b ) 4 } , which strongly 


suppresses the trapped particle contribution to w, when I co" \ > > w B , comes 
from the factor exp {- V| 2 /V T 2 } in (19a). 

To estimate the maximum growth rate, one sets w r % - u p and maximizes 
(23) with respect to «" . For w" = 2" % <^ B , 

- L 3 s«p(V ^) 2 


It is interesting to study the dependence of ( ^ on the amplitude of the 
whistler. If w T 2 is approximately proportional to the width of the trapping 
region in velocity space, then ~ and ( Wj) ~ e"' 4 . The growth rate of 
the instability therefore decreases with increasing e = Bj/B 0 . This result is 
not valid for extremely small e , because of the assumption that w. < < . 

Finally, consider the equilibrium established when a finite amplitude 
Whistler is "turned on" in an initially Maxwellian magnetoplasma and then 
allowed to evolve. The approximate asymptotic distribution function for trapped 
particles obtained 2 by the present authors for such a process is 


•x. 



'X 

*x 




~ 2^0 ^0 



This distribution shows no variation with W, and 3f T /Bv is less than zero for 
all v. In this case trapped particles help stabilize the plasma against perturba- 
tions having I w "l ~ w B . As in the preceeding case, the trapped particle con- 
tribution to the imaginary part of w becomes negligible when I I >> <^ B . 
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APPENDIX A 


We consider the effect of the transformation k 1 k + Alt, '(1<) (X )" (k + Alt) 

= - w" (It) on the quantity < g> k , when untrapped particles are ignored, Here g is 
any single valued function of v and YV only. It follows from (22) and (lib) that <o" 
occurs only in the combination w"/o> B in < g, ;i , so that <g> k - G (It, «" /<w B ). 

When | w"l is comparable to " B , Alt % -2 w (d w" /dlt)' 1 ~ O (e M ), Hence 
G (It + Alt, ~fo"/a> B ) ~ G (It, - <w"/to B ), to lowest order in e,. Using this fact, and 
again referring to (22), one finds that the sole effect of the transformation 
w" - - w " is to replace n by -n in each of the terms which contribute to the 
summation in (lib). Thus, to show that the trapped particle contribution to <g > k 
is an even function of o>", it is only necessary to show that |Aj ‘ = |A_J 2 . 


Equation (12c) may be written in the form 


S(v, w, t + , t) % y 


q ( m , + /r) ( _1 ) ,r,M cos 


t + -a 
(2m (■ 1) I ?.n — — + - 


(1 + q 2m+ 1) (2m + 1) 


m= o 


The function exp [ i k S] is periodic in [ 2 it (t - t + )/T + // /2 J , and may be ex- 
panded in a Fourier series. The coefficients of this expansion are, after using 
(6b), 

B n = e * in,,/2 A n 

But B n B_ n , because S is an even function of [2 n (t - t + )/T + n/2 \ , hence 

A„ a e ,nw A.„ 
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APPENDIX B 


A straightforward change of variables in (19a), using (1), (13), and the 
definition 

£ - + k o z + 'Pj > 

yields 





n eT k O H ( v2 " n o 2/k 0 2 ) 

exp 

2ft 0 n e 77(v)V T 2 


(v 2 ~ OW) 


S(l/r - 0) 5(0 


(Bl) 


f T must be a function of v, W or v, r only, i.e., the density of particles in phase 
space must be uniform along the particle orbits. The distribution described by 
(Bl) satisfies this requirement, even though it appears to have a third independ- 
ent variable. In fact, the orbit with 1/r = 0 consists of a single point on the 
line £ = 0. 

A more convenient representation of f T is obtained by first imagining the 
trapped particles to be evenly distributed aloi g a small but finite radius orbit 
near the bottom of the well, then allowing the radius of this orbit to go to zero. 
Thus 


'V. 



n eT ;: 0 H ( v2 ~ n 0 2/k 0 2 ) 
exp 

2n 0 n e 7 ?( V ) V T 


(v 2 ~ ft 0 2 / k 0 2 ) 

v T 2 


8 (1/r - 0) x C 


(B2) 


where C is determined by the condition 
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'+sin -1 (l/r) r+ s in -1 ( 1/r ) 

S (<5) (1/r 2 - sin 2 d£ = 2C (1/r 2 - sin 2 £)"“ d£ (B3) 

» e i n “ 1 / 1 •' r \ Cl n "* 1 / 1 ' r \ 


= 4 C K (1/r). 



C = r/ 4 K ( 1/r ) (B4) 

[See (22b)] . The factor of 2 on the right hand side of (B3) is due to the degen- 
eracy 8 of the coordinates r, v, After noting that K (1/r) ~ tt/2 for large r 
and inserting (B4) into (B2), one obtains the result quoted in (19b). 







■rs. 




/ 


REFERENCES 

1. R. F. Lutomirski and R. N. Sudan, Phys. Rev., 147, 156 (1966). 

2. P. J. Palmadesso and G. Schmidt, Phys. Fluids, To be Published June '971. 

3. C. S. Roberts and S. J. Buchsbaum, Phys. Rev., 135A , 381 (1964). 

4. M. J. Laird and F. B. Knox, Phys. Fluids, 8:, 755 (1965). 

5. W. L. Kruer, J. M. Dawson, and R. N. Sudan, Phys. Rev. Letters, 23, 838 

(1969). 

6. M. V. Goldman, Phys. Fluids, 13, 1281, (1970). 

7. B. B, Kadomtsev, A. B. Mikhailovskii, and A. V. Timofeev, Soviet Physics, 
JETP (English Translation), 1577 (1965). 

8. It should be noted that the coordinates r, v, and £ are degenerate in the 
sense that each point in r, v, g space represents two points in the original 
phase space; e.g., the two small circles in Figure 1 mark points with the 
same r, v, 


23 




•• 





have r > 1; untrapped particles haver < 1. cr - Sgn(V|| + figAg). 



